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Abstract 



The original annual-modulation region, singled out by the DAMA/Nal 
experiment for direct detection of WIMPs, is extended by taking into account 
the uncertainties in the galactic astrophysical velocities. Also the effect due to 
a possible bulk rotation for the dark matter halo is considered. We find that 
the range for the WIMP mass becomes 30 GeV < m v < 130 GeV at \-a C.L. 
with a further extension in the upper bound, when a possible bulk rotation of 
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the dark matter halo is taken into account. We show that the DAMA results, 
when interpreted in the framework of the Minimal Super symmetric extension 
of the Standard Model, are consistent with a relic neutralino as a dominant 
component of cold dark matter (on the average in our universe and in our 
galactic halo). It is also discussed the discovery potential for the relevant 
supersymmetric configurations at accelerators of present generation. 
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Typeset using REVTeX 



I. INTRODUCTION 



The DAMA/Nal Collaboration has recently reported the indication of an annual- 
modulation effect in a direct search experiment for WIMPs JT|-§[. In Ref. || it has been 
shown that a statistical (maximum-likelihood) analysis of the data concerning to a total 
exposure of 19,511 kg x day provides, at a 2— a C.L., a well delimited region in the plane 
Px a scaiar e0n ' ) ~~ m x> where m x is the WIMP mass, p x is its local (solar neighbourhood) density 
and Cgcata is the WIMP-nucleon scalar elastic cross section. Obviously, the location of 
the annual-modulation region depends on the functional form adopted for the speed distri- 
bution of the dark matter particles, and on the values assigned to the galactic astrophysical 
velocities. In Refs. |l]-[3| a standard Maxwellian distribution was taken for the WIMP speed 
distribution with a root mean square velocity t> rms = 270 Km s -1 and a WIMP escape 
velocity in the halo v esc = 650 Km s -1 . 

The DAMA/Nal results on annual-modulation were analyzed in Refs. f4j-|| in terms of 
relic neutralinos. It was proved that the DAMA/Nal data are compatible with a neutralino 
as a major component of dark matter in the universe. In Refs. it was also shown 

that a significant number of the relevant supersymmetric configurations are compatible with 
supergravity schemes ||, are explorable at accelerators and/or by indirect searches for relic 
particles (up-going muons at neutrino telescopes and antiprotons in space). 

In the present paper we discuss: i) the extension of the DAMA/Nal annual-modulation 
region in the plane p x o r s ^ 1 f r con ' ) - m x , when the uncertainties of the galactic astrophysical 
velocities are taken into account, and ii) the ensuing consequences of this extension in terms 
of properties for relic neutralinos and for the related supersymmetric configurations ||10|| . 

II. EXTENDED ANNUAL MODULATION REGION 

The differential rate for WIMP direct detection is given by 

dR p x r v ^ da , „ . 

- Nt ^2L / d #f(fl v —{v, E R ) , (1) 



dE R w dE R 

where is the number of the target nuclei per unit of mass, v and f(v) denote the WIMP 
velocity and velocity distribution function in the Earth frame (v = \v\) and dajdE R is 
the WIMP-nucleus differential cross section. The nuclear recoil energy is given by Er = 
m^ ed v 2 (l — cos 9*)/ mx, where 9* is the scattering angle in the WIMP-nucleus center-of-mass 
frame, is the nuclear mass and m re d is the WIMP-nucleus reduced mass. The velocity 
v m { n is defined as t> min = {ra^E R / 2ra^ e ^)^ and t> max is the maximal (escape) WIMP velocity 
in the Earth reference frame. Eq.(|TJ) refers to a mono-atomic detector; the generalization to 
our case of Nal is straightforward. 

In principle, the differential WIMP-nucleus cross section is composed of a coherent part 
and a spin-dependent one. Making the safe assumption that the former one is dominant 
over the latter, we may write 

^ ^ F\ Q ), (2) 



dE R E^ 
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where o"o is the point-like scalar WIMP-nucleus cross section, E™ x is the maximum value 
of Er and F(q) denotes the nuclear form factor, expressed as a function of the momentum 



transfer q 



WIMP-nucleon cross section v^aiac 



2mjyEji. o"o may be conveniently rewritten in terms of the scalar 



Co 



( l+m x /m p \ 2 2 

I ~, I ■* i U, 

\l + m x /m N J 



(nucleon) 
scalar ' 



(3) 



where m p is the proton mass and A is the mass number of the nucleus. 

For the nuclear form factor in Eq.(|2]) we use the Helm parameterization of the scalar 
form factor |TT1JT^] : 



F(q) 



oJlOo) 

3 exp 

qr 



1 



2 2 

s q 



(4) 



where s ~ 1 fm is the thickness parameter for the nucleus surface, r = (r 2 — 5s 2 ) 1 / 2 , 
r = 1.2 A 1 / 3 fm and ji(qro) is the spherical Bessel function of index 1. 

Once a functional form for f(v) is chosen and specific values are assigned to the relevant 
astrophysical velocities, the quantity Px a scaiar e0n ' ) ma y be extracted from measurements of 
the differential rate dR/dEn function of m v . 

Here, for f(v) we take the standard Maxwell-Boltzmann distribution (as in an isothermal 
spherical model for the halo) with a finite escape velocity, i.e. in the galactic rest frame we 
write 



3/2 



,gah 



N 



2irv 2 



exp 



"rms / 



where the normalization factor is 

N 



i -i 



erf(z) =z exp(—z 2 

7T 



(5) 



(6) 



is related to the asymptotic 
The measured rotational 



with z 2 = 3f 2 sc /(2t> 2 ms ). In the isothermal spherical model, t> rms 
value v oo of the rotational velocities by the relation t> rms = yfi 
velocity remains almost flat (to roughly 15%) between 4 Kpc and 18 Kpc; here, we identify 
Voo with the rotational velocity of the Local System at the position of the Solar System 
v (r Q ) = Vq, whose physical range will be shortly discussed. 

To make use of Eq. ([]]), Eq.(|5|) has to be transformed to the rest frame of the Earth, 
which moves through the Galaxy with a velocity 



v (B = v Q + v orh cos 7 cos[uj(t - t )], 
in the azimuthal direction. v & is given by 

Vq = v o + 12 km s _1 . 



(7) 



(8) 



In these expressions the speed of 12 km s 1 stands for the motion of the Solar System with 
respect to the Local System, f or b = 30 km s -1 denotes the Earth orbital speed around the 
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Sun, the angle 7 ~ 60 is the inclination of the Earth orbital plane with respect to the 
galactic plane and u = 27r/365days, t = June 2 nd [13,14 . 



In Refs. ||T]-§J the DAMA/Nal annual-modulation region was extracted from the experi- 
mental data, by using a maximum likelihood method and by taking for the velocities vq and 
v esc the following values: vq = 220 km s _1 and f esc = 650 km s" 1 (the relevant region is one 
of the domains displayed in Fig. 1). The ensuing implications for relic neutralinos, derived 
in j§-§] , refer to the same set of astrophysical velocities. 

In the present paper we extend the analysis of the annual-modulation region, by consid- 
ering the physical ranges associated to v and t> esc [|T3|— P~7|] : 

v = (220 ± 50) km s" 1 (90% C.L.) , (9) 



Vesc = (450 -=- 650) km s" 1 (90% C.L.) . (10) 

The statistical method for the extraction of the annual-modulation region is the same 
as the one employed in Refs. [0-0], with a lower bound m x > 25 GeV. 

On the basis of the analytical properties of the time-modulated part of the detection rate 



14], one expects that a variation in vq induces a sizeable modification in the range of m x , 
without affecting p x c^^° significantly. Indeed, an increase (decrease) in v is expected to 
extend the original DAMA/Nal region toward lower (larger) WIMP masses for kinematical 
reasons; Px^scaiar cannot appreciably change, since it acts as a normalization factor which 
is determined by the size of the average detection rate. Also, no sizeable variation of both 
m x and p^sadar*"^ * s expected from a variation in t> esc , since the escape velocity provides 
a cut-off in the integral of Eq.(0) applying only on the flat tail of the velocity distribution 
/('•)• 

We report in Fig. 1 the (2-er C.L.) annual-modulation regions extracted from the 
DAMA/Nal data of Ref. || (total exposure of 19,511 kg x day) for sets of values for 
Vq and v csc , which bracket the ranges given in Eqs.(§ - [10]). The 1-er ranges for the relevant 
quantities are given in Table 1. Notice that in Fig. 1 the density p x is given in units of 0.3 
GeV cm" 3 , i.e. p° x 3 = pj (0.3 GeV cm" 3 ). 

From the features displayed in Fig. 1 we notice that, as anticipated above: i) the location 
of the annual-modulation region is rather sensitive to the velocity Vq, with a sizeable exten- 
sion in the range of m x , but with small variations in Px a scaiar e ° n ' ) j ii) the annual-modulation 
region is essentially independent of v esc . From the results of Table I we conclude that the 
uncertainties in Vq extend the range of m x , as singled out by the DAMA/Nal annual- 
modulation data, to (at 1-er C.L.) 

30 GeV < m y < 130 GeV . (11) 

Taking into account the whole ranges of v and v esc (see Eqs. (^fLOf)), we find the annual- 
modulation region which is depicted in Fig. 2. This region, which envelops the domains 
displayed in Fig. 1, will be hereafter referred to as region R m . 

It is worth noticing that this annual-modulation region might be further extended, if a 
possible bulk rotation of the dark matter halo is introduced. Unfortunately, halo models 
which take realistically into account this phenomenon are not yet available; however, effects 
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of rotation of the dark matter halo on direct detection rates and ensuing upper bounds on 
cross sections have been addressed in some extreme models [|l8|,|rjj . To obtain an estimate of 
the effects of a possible rotation of the isothermal sphere one can consider a class of models 
which describe the fastest rotating steady state by means of the following recipe: 

UiV )" vT x <0 (12) 



/> sal ) 




vT l > o 



,gal 



< 



(13) 



where vT l is the azimuthal component of v. 

Then one can combine the above functions in order to deal with a more general family 
of distribution functions f T ot(v gal ), defined as fT8fl : 



/rot(^ gal ) = af + (v^) + (1 - a)/_K al ) 



,gal 



(14) 



where a is related to the dimensionless galactic angular momentum parameter A by the 
relation A = 0.36|a — 0.51. In order to be consistent with the available extensive numerical 



work on galaxy formation, A should not exceed the value 0.05 [21 



Adopting the expression in Eq. (jT4|) for the velocity distribution function, one finds that 
the annual-modulation region may be extended to the domain depicted in Fig. 3, and the 
relevant range of m x would become (at 1-er C.L.) 



30 GeV < m y < 180 GeV. 



(15) 



Although the possible occurrence of a bulk rotation of the dark matter halo is a quite 
interesting possibility deserving further investigation, all subsequent analyses of the present 
paper will refer to the region R m , given in Fig. 2. 



III. SUPERSYMMETRIC MODEL 

The WIMP candidate considered in this paper is the neutralino, defined as the lowest- 
mass linear superposition of photino (7), zino (Z) and the two higgsino states (H°, H%) 

X = ai7 + a 2 Z + a 3 H° + a A H^ . (16) 

To classify the nature of the neutralino it is useful to define a parameter P = a\ + a\\ 
hereafter the neutralino is called a gaugino, when P > 0.9, is called mixed when 0.1 < P < 
0.9 and a higgsino when P < 0.1. 

The theoretical framework adopted here is the Minimal Supersymmetric extension of 
the Standard Model (MSSM) |p2[ , which conveniently describes the supersymmetric phe- 
nomenology at the electroweak scale, without being constrained by too strong theoretical 
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assumptions. Specific details of the scheme employed in this paper are given in Ref . || ; here 
we only recall a few essentials. 

The large number of free parameters inherent in the model is reduced to six independent 
ones, by imposing a few assumptions at the electroweak scale: i) all trilinear parameters 
are set to zero except those of the third family, which are unified to a common value A; ii) 
all squarks and sleptons soft-mass parameters are taken as degenerate: mj, = = mo, 
iii) the gaugino masses are assumed to unify at Mgut, and this implies that the U(l) and 
SU(2) gaugino masses are related at the electroweak scale by Mi = (5/3) tan 2 Q W M 2 . 

The six independent parameters are taken to be: M 2 , /i, tan f3, m^, m , A, where p is the 
Higgs-mixing parameter, tan (3 is the ratio of the two Higgs vacuum expectation values, and 
rriA is the mass of the neutral pseudoscalar Higgs boson. To get the scatter plots shown in 
this paper, the supersymmetric space has been scanned, by varying the parameters in the 
following ranges: 10 GeV < M 2 < 1000 GeV, 10 GeV < < 1000 GeV, 80 GeV < m A < 
1 TeV, 100 GeV < m < 1 TeV, -3 < A < +3, 1 < tan/3 < 50. 

Our supersymmetric parameter space is further constrained by all the experimental limits 
obtained from accelerators on supersymmetric and Higgs searches. The latest data from 
LEP2 on Higgs, neutralino, chargino and sfermion masses are used |24[]. Also the constraints 
due to the b — > s + 7 process [^5|,^J have been taken into account (see Ref. M for the 
theoretical details). 

The supersymmetric parameter space has also been constrained by the requirement that 
the neutralino is the Lightest Supersymmetric Particle (LSP), i.e. regions where the gluino 
or squarks or sleptons are lighter than the neutralino have been excluded. 

One further constraint is due to the requirement that the neutralino relic abundance does 
not exceed the cosmological bound, derivable from measurements of the age of the Universe 
27[ and of the Hubble constant We have adopted here a conservative upper bound, 



Vt x h < 0.7 (h is the usual Hubble parameter, defined in terms of the present-day value Hq 
of the Hubble constant as h = Hq/(100 km s -1 Mpc -1 )). 



The neutralino relic abundance is calculated here as in Ref. |29| . Inclusion of coannihila- 
tion effects |30 in the calculation of Q x h 2 are not necessary here, since the instances under 
which these effects might be sizeable are marginal for the supersymmetric configurations 
concerning the DAMA/Nal data. 

Whenever we have to evaluate the neutralino local density, we employ the rescaling 
procedure. This rescaling consists in assuming that the neutralino local density p x may 
be taken as p x = t^pi (pi is the total local density of non-baryonic dark matter), with 
£ = min (1, f2 x /i 2 /(fi/i 2 ) min ), i.e., p x may be set equal to pi only when £l x h 2 is larger 
than a minimal value (Qh 2 ) min , compatible with observational data and with large-scale 
structure calculations; otherwise, when fl x h 2 turns out to be less than (flh 2 ) m i n , and then the 
neutralino may only provide a fractional contribution Q x h 2 /(Qh 2 ) m m to Qh 2 , p x is reduced 
by the same fraction £ = Q x h 2 /(Qh 2 ) min as compared to p\. The value to be assigned 
to (Qh 2 ) min is somewhat arbitrary, in the range 0.01 < (Qh 2 ) min < 0.2. In the present 
paper, whenever we have to apply rescaling, we use the value (Qh 2 ) min = 0.01, which is 
conservatively derived from the estimate f2 ga iactic ~ 0.03. 

In Fig. 2 we display our results for p° ,s cr s ^ 1 ^ r eon ' ) in the form of the scatter plot which is 
derived by scanning the supersymmetric parameter space over the grid defined above. For 
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the evaluation of cr^^J we use the expressions given in Ref. 0. We notice that a host of 
configurations fall inside the region R m ; thus, the annual-modulation region is compatible 
with supersymmetric configurations currently allowed by accelerator constraints. 

IV. INTERPRETATION OF THE ANNUAL MODULATION DATA IN TERMS 

OF RELIC NEUTRALINOS 

Let us turn now to the implications of the DAMA/Nal experimental data, once these 
are interpreted in terms of relic neutralinos. The cosmological properties are examined first; 
other prominent features of the relevant supersymmetric configurations will be examined 
afterwards. 



A. Neutralino cosmological properties 

For the derivation of the relic neutralino properties as regards its local density p x as 
well as its contribution to the average relic abundance Q x h 2 compatible with region R m , we 
adopt the straightforward procedure, outlined in Ref. ||, which does not require any use of 
rescaling in the local neutralino density. The method is the following: 

1) We evaluate 0^1^°^ an d £l x h 2 by varying the supersymmetric parameters over the 
grid and with the constraints defined in the previous section. 

2) For any value of [Px cr scaiar e ° n ' ) ]expt compatible with the region R m we calculate p x as 

given by p x = [Px a scaiar 0On ^]cxpt/^criar eOn ' ) an d restrict the values of m x to stay inside the region 
Rm displayed in Fig. 2. 

3) The results are then displayed in a scatter plot in the plane p x vs. Q x h 2 . 
Examples of our results are given in Figs. 4 a-c for a few experimentally allowed values 

of 

Px^icaiar ^> which bracket the range implied by the region R m of Fig. 2. Sects, a— c of 
Figs. 4 refer to the values [p^' 3 ^scaiar e ° n ' ) ]cxpt = 4 • 1CT 9 nbarn, 6 • 10~ 9 nbarn and 8 • 1CT 9 
nbarn, respectively. 

The two horizontal lines delimit the physical range 0.1 GeV cm -3 < p\ < 0.7 GeV cm -3 
for the total local density of non-baryonic dark matter. This (rather generous) range has 



been established by taking into account a possible flattening of the dark matter halo [31,32 



and a possibly sizeable baryonic contribution to the galactic dark matter [33]. The solid 
vertical lines delimit the cosmologically interesting the range 0.01 < Vt x h 2 < 0.7. The two 
vertical dashed lines delimit the range: 0.02 < Q x h 2 < 0.2, which represents the most 



appealing interval. Indeed, recent observations and analyses [|4[] favour 0.1 < ^matter < 0.4. 



Combining this range with the one for h: 0.55 < h < 0.80 and requiring that a cold 
dark matter candidate (such as the neutralino) supplies ~ (80-90)% of fi ma tter, we actually 
obtain: 0.02 < ficDM^ 2 < 0.2. 

The two slant dot-dashed lines delimit the band where linear rescaling procedure for the 
local density is usually applied. In Figs. 4 a-c the upper dot-dashed line would refer to a 
rescaling with (Qh 2 ) min = 0.01, the lower one to the value (Qh 2 ) min = 0.2. However, notice 
that in the derivation of the scatter plot of Figs. 4 a ~ c > no use of rescaling for p x is made. 
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With the aid of this kind of plot we can classify the supersymmetric configurations 
belonging to region R m into various categories. Configurations whose representative points 
fall above the maximum value p x = 0.7 GeV cm~ 3 have to be excluded (we remind that 
those providing an Vt x h 2 > 0.7 are already disregarded from the very beginning). Among 
the allowed configurations, those falling in the region inside both the horizontal and solid 
vertical lines (called A hereafter) are very appealing, since they would represent situations 
where the neutralino could have the role of a dominant cold dark matter component; even 
more so, if the representative points fall in the subregion (called B hereafter) inside the 
vertical band delimited by dashed lines. Configurations which fall inside the band delimited 
by the slant dot-dashed lines denote situations where the neutralino can only provide a 
fraction of the cold dark matter both at the level of local density and at the level of the 
average fl. Configurations above the upper dot-dashed line and below the upper horizontal 
solid line would imply a stronger clustering of neutralinos in our halo as compared to their 
average distribution in the Universe. 

It is worth noticing a few important properties of the scatter plots shown in Figs. 4 a-c: 

1) The scatter plots display a correlation between p x and Q x h 2 . This feature is expected 
on the basis of the following properties: i) Q x h 2 is roughly inversely proportional 
to the neutralino pair annihilation cross section, ii) at fixed [p x ^scaii« ]expt> P x * s 
inversely proportional to (T^^° , iii) the annihilation cross section and cr^\^° are 
usually correlated functions (i.e., they are both increasing or decreasing functions of 
the supersymmetric parameters). 

2) The domains covered by the supersymmetric configurations in regions A and B are 
slightly larger for smaller values of [Px^scaiar ]expt- This feature follows from the fact 
that crg^ 1 r eon ' 1 is bounded from above by accelerator limits (mainly because of lower 
bounds on Higgs masses); this implies for p x a lower bound, which however is less 
stringent at lower values of [Px a scater° n ^expt- 

We can conclude that the DAMA/Nal data are compatible with a relic neutralino as a 
major component of dark matter in the universe. 



B. Other properties of the supersymmetric configurations of set S m 

Apart from the cosmological properties previously discussed, one of the most interesting 
questions is whether the supersymmetric configurations, whose representative points fall 
inside the region R m of Fig. 2 (hereafter denoted as configurations of set S m ), are explorable 
at accelerators of the present generation. This point was investigated in Ref. M in the case 
of the original DAMA/Nal annual-modulation region. Here we extend our considerations 
to the larger set S m derived from region R m . 

The supersymmetric configurations of set S m are constrained by the inequalities 

r (nucleon) -I r (nucleon) -i 

[Px^scalar Jcxpt , (nucleon) , [Px^scalar Jexpt (-\*7\ 

- ^ °scalar — ~ \ L 1 ) 

Px> max rXi m i n 
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where p Ximin = £p«,min = £• 0.1 GeV cm 3 and p x , max = £,pi, max = £• 0.7 GeV cm 3 , and 
[Px^sc^s^ lexpt is an Y value inside region R m as a function of m x . 

Eq.(|T7|) implies interesting correlations among the supersymmetric parameters, as dis- 
played in Figs. 5-7. The first of these figures shows the scatter plot of tan (3 versus the mass 
rrih of the lightest neutral CP-even Higgs boson (we recall that rrih may be derived from the 
supersymmetric parameters listed in Sect. III). The correlation displayed in Fig. 5 between 
tan (3 and rrih is implied by the fact that &s™te,r° is usually dominated by Higgs-exchange 
amplitudes, and these are in turn large for large values of tan/3 and small values of rrih. 
As is apparent in Fig. 5 a number of supersymmetric configurations of set S m could still 
be explorable at LEP2, but the others require investigation at a high luminosity Fermilab 
Tevatron |35| , |36| or at LHC Fig. 6 shows a correlation between tan/3 and the mass of the 
lightest top-squark m^. The reasons for this feature are more involved, and are discussed 
in Ref. §. 

A last plot, providing tan/? versus m x , is given in Fig. 7. Again, a correlation shows up 
here for not too large values of tan (3. LEP2 can only provide an investigation up to m x ~ 
50 GeV. The exploration up to ~ 125 GeV can be performed with Tevatron upgrades (under 
favorable conditions), while for higher values of m x LHC is needed. 



V. CONCLUSIONS 

We have investigated how the original DAMA/Nal annual-modulation region in the plane 
Px cr scaiar e0n ' ) ~~ m x * s extended when the uncertainties in the galactic astrophysical velocities 
are taken into account. One of the most noticeable consequences is that the range for the 
neutralino mass becomes 30 GeV < m Y < 130 GeV at 1-er C.L., with a further extension 
in the upper bound, when a possible bulk rotation of the dark matter halo is taken into 
account. However, the extension of the modulation region implies only slight modifications 
in the correlations among the supersymmetric parameters in the MSSM scheme. 

We have found that a number of supersymmetric configurations singled out by the 
DAMA/Nal results have cosmological properties compatible with a relic neutralino as a 
dominant component of cold dark matter (on the average in our universe and in our galactic 
halo). It has also been discussed the discovery potential for the relevant supersymmetric 
configurations at accelerators of present generation. 
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Note added 

During the completion of this work there appeared two preprints: [hep-ph / 9903467] (by 
L. Roskowsky) and [hep-ph / 9903468| (by M. Brhlik and L. Roszkowski), where some consid- 
erations on the effect of astrophysical velocities uncertainty on WIMP direct searches are 
presented. 
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TABLES 



TABLE I. Results of the maximum likelihood method when varying the velocity parameters 
(vo and w csc ) of the WIMP velocity distribution. The values of m x and v^^^^ correspond to the 
minima of the maximum likelihood method (with 1— a error bars). 
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FIG. 1. Annual-modulation regions singled out by the DAMA/Nal experiment [3] in 
the plane Py 3 Cscaiar° n ^ ~ m x f° r var i° us values of velocities v and v e 



v esc = 650 km s (solid line); Vo = 220 km s , v esc 
vq = 170 km s -1 , v esc = 550 km s _1 (dotted line); vq 
(dashed line). The countour line for v = 220 km s -1 , v es 



v = 220 km s -1 , 
= 450 km s _1 (dot-dashed line); 
= 270 km s -1 , v esc = 550 km s _1 
= 550 km s~ x is undistinguishable 



on the plot from the solid line. The quantity p° 3 is the local neutralino matter density in 
units of 0.3 GeV cm" 3 , i.e. p°' 3 /(0.3 GeV cm" 3 ). 
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FIG. 2. The contour line delimits the annual-modulation region obtained by varying the 
velocities vq and v esc over the ranges given in Eqs. (|9| - PHI). The scatter plot represents 
the theoretical predictions of a generic MSSM, as described in Sect. III. Different symbols 
identify different neutralino compositions: circles stand for a higgsino, crosses for a gaugino 
and dots for a mixed neutralino. 
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FIG. 3. Annual-modulation regions when a possible bulk rotation of the dark matter 
halo is taken into account. The contour lines refer to the following halo parameters: Vq = 220 
km s -1 , non rotating halo (solid line); vq = 220 km s" 1 , co-rotating halo with spin parameter 
A = 0.05 (dotted line); v = 220 km s -1 , counter-rotating halo with A = 0.05 (dashed line); 
Vq = 170 km s -1 , co-rotating halo with A = 0.05 (dot-dashed line); v = 270 km s -1 , 
counter-rotating halo with A = 0.05 (long-dashed line); The escape velocity is fixed at the 
value v esc = 550 km s" 1 for all the contours. 
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FIG. 4a. The neutralino local density p x , calculated for fixed values of [p x cr^^ n ] exp t, 
is plotted versus the neutralino relic abundance Q x h 2 . For any value of [p x ^&scaiar n }expt the 
neutralino mass is restricted to the range dictated by the contour line of Fig. 2. The different 
panels of the figure correspond to the following choices: (a) [p°'V" c u a f a c r on ] eX p t = 4 • 1CT 9 nbarn 
and 28 GeV < m x < 88 GeV; (b) [p x 3 cr™£™%xpt= 6 • 1(T 9 nbarn and 25 GeV < m x < 131 
GeV; (c) [p° x (y^° n ]e Xp t= 8 • 10~ 9 nbarn and 25 GeV < m x < 156 GeV. The two horizontal 
lines delimit the physical range for the local density for non-baryonic dark matter. The two 
solid vertical lines delimit the interval of Q x h 2 of cosmological interest. The two vertical 
dashed lines delimit the preferred band for cold dark matter. The two slant dot-dashed 
lines delimit the band where linear rescling procedure is ususally applied. 
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FIG. 4b. The neutralino local density p x , calculated for fixed values of [p x ^&fcai&° n ]expt, 
is plotted versus the neutralino relic abundance Q x h 2 . For any value of [p x 3 o f^iar m ]expt the 
neutralino mass is restricted to the range dictated by the contour line of Fig. 2. The different 
panels of the figure correspond to the following choices: (a) [p°' 3 ^ c u a f a e ° n ]e:rpt= 4 • 1CT 9 nbarn 
and 28 GeV < m x < 88 GeV; (b) [p x 3 cr™£™%xpt= 6 • 10" 9 nbarn and 25 GeV < m x < 131 
GeV; (c) [p° x 2 l <" a f£ on W= 8 • 10 -9 nbarn and 25 GeV < m x < 156 GeV. The two horizontal 
lines delimit the physical range for the local density for non-baryonic dark matter. The two 
solid vertical lines delimit the interval of Q x h 2 of cosmological interest. The two vertical 
dashed lines delimit the preferred band for cold dark matter. The two slant dot-dashed 
lines delimit the band where linear rescling procedure is ususally applied. 
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FIG. 4c. The neutralino local density p x , calculated for fixed values of [p x cr^^ n ] exp t, 



x 

is plotted versus the neutralino relic abundance VL x h 2 . For any value of [p x 3 o '^iar° n ] expt the 
neutralino mass is restricted to the range dictated by the contour line of Fig. 2. The different 
panels of the figure correspond to the following choices: (a) [p°' 3 ^ c u a f a e ° n ]e:rpt= 4 • 1CT 9 nbarn 
and 28 GeV < m x < 88 GeV; (b) [p x 3 cr™£™%xpt= 6 • 10" 9 nbarn and 25 GeV < m x < 131 
GeV; (c) [p° x 2 '(r^° n } expt = 8 • 10" 9 nbarn and 25 GeV < m x < 156 GeV. The two horizontal 
lines delimit the physical range for the local density for non-baryonic dark matter. The two 
solid vertical lines delimit the interval of Q x h 2 of cosmological interest. The two vertical 
dashed lines delimit the preferred band for cold dark matter. The two slant dot-dashed 
lines delimit the band where linear rescling procedure is ususally applied. 
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FIG. 5. Scatter plot for set S in the plane - tan/3. The hatched region on the right 
is excluded by theory. The hatched region on the left is excluded by present LEP data at 
y/s = 189 GeV. The dotted and the dashed curves denote the reach of LEP2 at energies y/s 
= 192 GeV and y/s = 200 GeV, respectively. The solid line represents the 95% C.L. bound 
reachable at LEP2, in case of non discovery of a neutral Higgs boson. 



20 




FIG. 6. Scatter plot for set S in the plane - tan (3. The hatched region is excluded 
by LEP data. 
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FIG. 7. Scatter plot for set S in the plane m x - tan (3. The hatched region on the left is 
excluded by present LEP data. The dashed and the solid vertical lines denote the reach of 
LEP2 and TeV33, respectively. 
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